PROTEUS MIRABILIS SWARM-COLONY DEVELOPMENT WITH DRIFT 



PHILIPPE LAURENgOT AND CHRISTOPH WALKER 



Abstract. We prove a global existence result for a model describing the swarming phenomenon 
of the bacterium Proteus mirabilis. The model consists of an ordinary differential equation cou- 
pled with an age-structured equation involving nonlinear degenerate diffusion and an additional 
drift term. 



1. Introduction 

Proteus mirabilis is a bacterium that is widely distributed in soil and water in the natural 
environment. It is also found in the intestinal tract of many mammals, including human. Broth 
cultures of Proteus mirabilis consist of small swimmer cells, but produce a morphologically and 
physiologically distinct cell type, called swarmer cells, when inoculated on a solid surface. This 
process is referred to as "differentiation" and is crucial for the pathogenesis of these bacteria during 
urinary tract infections caused, e.g., by long-term urinary catheterization. While swimmer cells go 
through a prototypical cell division process and are immobile, swarmer cells age and increase in size. 
Swarmer cells can group together to build multicellular "rafts" that, when of sufficient biomass, 
are capable of translocation. This leads to a migration phase during which swarmer cells may 
also dedifferentiate again into swimmer cells. Once the biomass falls below the critical threshold, 
movement ceases initiating a consolidation phase. This oscillation between phases of motion of 
swarmer cells and consolidation to the swimmer state leads to an interesting bull's-eye-patterned 
biofilm. 

Different mathematical models were proposed in order to describe the swarming of Proteus 
mirabilis. In this article we focus on a model that was presented and numerically analyzed in 
[2] and later, in slightly modified form, in [31 [TC]. The model involves the swimmer cell density 
V = v{t, x) in dependence of time t and spatial position x and the swarmer cell density u = u{t, a, x), 
where the variable a models cell age. The equations under consideration are 

dtu + daU ^ divx{D{A)\/xU + uE{A,v)\/xA) - fi{a)u , (t, a, e (0, oo)^ x , (1.1) 



dtv = [g{v) - i{v)) V + b{a) ^(a) u{t, a, x) da , {t, x) e (0, oo) x f2 , (1.2) 

where 

/•OO 

A(t, x) / \{a)u{t,a,x)da , (t, x) £ (0, cx)) x f2 , (1.3) 

^0 

subject to the boundary conditions 

u{t,Q,x) ^ S,{v{t,x))v{t,x) , (t,x) e (0,oo) x 17 , (1.4) 

D{K)dyU + u E{A, v) d^A ^ , {t, a, x) e (0, oo)^ x dn , (1.5) 

and the initial conditions 

m(0, a, x) = u°(a, x) , v{0, x) = v''\x) , (a, x) £ (0, oo) x 17 . (1.6) 
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Here, C K.^ is an open and bounded set with smooth boundary dH, and i' — (ly^ , . . . , z^^) denotes 
the outward normal unit vector field to dfl. The function A given by Ijl.Sp represents the total 
motile swarmer cell biomass. Since increase in size of swarmer cells is exponential with increase in 
age, the function A appearing in the definition of A is often taken in the form 

A(a) = mol[ao,oo)(a)e°^'' , 

where ag > is the minimal age of swarmer cells required to participate actively in group migration. 
The parameter r is the average time it takes a cell to subdivide, and mo > is a constant. 

Equation (|l.ip expresses the change in time of swarmer cells of a given age a. The diffusivity 
D depends on A and is zero for A = or, more generally, for A small. The explicit appearance of 
the "drift" term divx(uE{A,v)'Vx^) on the right-hand side of Ijl.ip is a novelty, although it has 
been already implicitly contemplated in the model of More precisely, in 3J the diffusion term 
is derived from isotropic random motion (instead of Fickian diffusion) and thus written in radially 
symmetric coordinates as {l/r)dr{rdr{DA{A)u)) which corresponds to the choice D — Da and 
E{A, v) = D'j^{A) in our case. As pointed out in only swarmer cells of a certain maturity can 
actively participate in group migration but nothing prevents young swarmers from being caught up 
in the flow and thus move with larger swarmers in the rafts. However, diffusion terms of the form 
diyx{D{A)Vxu) as considered in [HlIinilTB] reflect active movement of swarmers of any age, i.e. also 
of young swarmers. Therefore, we hypothesize that this term could be rather small (or even zero) 
and that migration of swarmers could be mainly due to the drift term div^, (uE{A, t;) Va;A) in which 
small (i.e. young) swarmers move but do not actively contribute to a raft's motility. Interesting 
would be, of course, to see whether numerical computations can support this hypothesis. The age 
dependent function fi in (|l.ip is the dedifferentiation modulus, which is higher for older swarmers 
than for younger ones. 

The change in time of the swimmer population is given by equation (|1.2p . The population grows 
exponentially, where often g{v) = in numerical simulations. Swimmer cells differentiate with 
rate ^{v) into swarmers of age leading to the age boundary condition (|1.4p . Usually, ^ is of the 
form ^(w) = ^{v)/t^ where ^(u) = for small and large values of respectively. The incorporation 
of a lag phase in swarmer cell production triggers the development of the consolidation phase after a 
swarm phase. It thus prevents a self-sustaining soliton caused by swarmers that dedifferentiate into 
swimmers immediately differentiating into new swarmers. This lag in the onset of differentiation 
was used in [3l |4l [16] . The integral term in (|1.2p represents dedifferentiation of swarmer cells into 
swimmer cells. In the numerical simulations in [31 HI [SI [1^] the function b is given by h{a) = e"/"^. 

As for further explanation of the model and for computational results regarding (|l.ip - (|1.6p we 
refer to [3l |4l [9l [Til [l^ a-nd the references therein. Existence results for (ll.ip - (|1.6p in the case of 
non-degenerate diffusion and E = can be found in '10*, while degenerate diffusion and E = 
was studied in [14J. 

The purpose of this article is to investigate mathematically equations (ll.ip - (|1.6p for degenerate 
diffusion with non-vanishing drift term. As already pointed out, the latter, in our opinion, could 
possibly be more important in the migration process than diffusion. It also generates additional 
difficulties in the mathematical analysis: indeed, while an Loo-estimate for u is readily obtained 
from p.ip when E = 0, such a bound does not seem to be available in the presence of the drift 
term (E ^ 0) and we only obtain a much weaker L In L-estimate on u, see Lemma 13.61 below. As a 
consequence of this lower regularity and the possible degeneracy of D and E, the diffusion and drift 
terms are not well-defined in (jl.ip and a weak formulation is required. The latter is introduced 
in Definition 12.11 and is somehow reminiscent of the definition of renormalized solutions for the 
Boltzmann equation [^ or parabolic equations (see, e.g., [2l[5l|6] and the references therein). 

Before stating precisely the assumptions on the data used in this paper and the results obtained, 
let us briefly outline the difficulties to be overcome and sketch our approach. First, the system 
(|l.ip - (|1.6p is of mixed type and features several nonlinearities. We will thus use a compactness 
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method, that is, first estabhsh the existence of solutions to a sequence of approximate problems 
and then pass to the limit as the approximation parameter converges to zero. Besides standard 
approximations (such as the positivity of D, the boundedness of E, and an additional linear 
diffusion in (jl.2p '). the approximation used herein relies upon the discretization of (jl.ip with 
respect to the age variable which leads us to a system of parabolic equations to which the abstract 
theory developed by Amann [r can be applied. Next, (|l.ip is a first-order transport equation with 
respect to the age variable and a degenerate parabolic equation with respect to the space variable, 
while (|1.2p features no spatial diffusion. The latter thus does not provide any smoothing effect 
which would guarantee the strong compactness for v needed to pass to the limit in the nonlinear 
terms E{A, v), g{v), and £,{v). Strong compactness for v can thus only result from that of the last 
term of the right-hand side of (|1.2p which requires the strong compactness for the age averages of 
u. Such a compactness property can only be deduced from (jl.ip but is hindered by the possible 
degeneracy of the diffusion term D{A)VxU for small values of A. Nevertheless, it holds true under 
the assumptions that A is bounded from below by a positive constant while D only vanishes when 
A vanishes. It is, however, the main obstacle to include the case D = in our analysis, see Remark 
14.61 As already mentioned, another difficulty to be faced is that the drift term diVrc{uE{A, u)V^A) 
only allows us to obtain an estimate of u in LlnL while E{A,v)W ^A belongs merely to L2, so 
that the drift term is not well-defined. Here, we take advantage of the Loo-boundedness of the 
age averages of u to set up a weak formulation which complies with the available regularity of u. 
Concerning compactness estimates for u, they rely on the just mentioned L In L-estimate derived 
from the specific structure of (|l.ip as well as a degenerate parabolic equation in the variables t 
and X satisfied by the age averages of u. While the former guarantees the weak compactness for u 
in Li, the latter provides the expected strong compactness on the age averages of u thanks to our 
assumptions on the data. 

We shall also remark that one can use different approaches to investigate the existence of solu- 
tions to age structured equations with (non-degenerate) diffusion, including integrated semigroups, 
perturbation arguments, or using solutions integrated along characteristics (e.g., see jTS l ITS l [20] 
and the references therein.) However, handling such equations by discretizing with respect to the 
age variable seems to be a novel approach. Let us also point out that the method used in [14j 
to tackle the case of age structure with quasi-linear non-degenerate diffusion is apparently not 
applicable in the present situation due to the additional drift term. 

Remark 1.1. In the particular case where A(a) = mo e"/"^, 6(a) = mi e"/"^, and ^{a) = m2 for 
some positive real numbers mo, mi, m2, and t, a closed system for the evolution of A and v can 
be derived from (|l.ip - (|1.2p and reads 

dtA = div, ( {D{A) + A E{A, v)) V, A) + - A , {t, x) e (0, 00) x , 

dtv = (g(„) _ ^(v))„ +!I^A , {t,x)^{0,oo)xn . 

mo 

The analysis performed in |16j actually focuses on this "reduced" system with the choice of 

DoA 



D{A) + AE{A,v) = 
where Dn, fc > 0. 



A + kv ' 



2. Existence 

Regarding the data in (|l.ip - (ll.6p we will assume that the following hypotheses hold: 

(hi) The functions $,5 G L^{R) n C^{R) are such that ^(s) = for s < and < ^(s) < g(s) 
for s > 0. 
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(/i2) The function b G C^([0, oo)) is non-decreasing with 6(0) = 1 and b{a) — > oo as a — > oo, and 
there exists a number Bq £ (0, oo) such that 

b{a + a) < {Boa + 1) b{a) , a > , a G (0, 1) . 

(/13) The function A G C^([0,oo)) is non-negative, satisfies £0 '■— inf A > 0, and there exists a 
number G (0, 00) such that 

A(a + a) < (Lo a + 1) A(a) , a > , a G (0, 1) , 

A(a - a) < {Lqu + 1) A(a) , a > a , a G (0, 1) . 

(/14) The function G Loo(0,oo) is non-negative, and there exists a number /3o G (l,oo) such 
that 

Ai(a) Ko) < Po A(a) < 6(a) , a > . 
(/15) The function D G C'^(R) is non-decreasing with D{r) > for r > and [r i— > {D{r)/r)^/^] G 
Li(0, 1). Moreover, for the function Ci, defined by Cl(f) := {D{r)/r)^^^ and Ci(0) = 0, we 
assume that f'/Ci G C([0, 00)) and put 

D'(r) 

Ki{R) sup — f-^- , i? > . 
0<r<i?, Ci('') 

(/ig) The function E G C'^(R^) is non- negative, and there is a function C2 G C^(R) such that 
(2(0) = 0, C2ir) > for r > 0, E^/C^ G C([0, 00) x [0, 00)), and 

^2iR) < Eir,s) < KsiR) C2W , ir,s) G [0,i?] x [0, i?] , 

for some constants K2{R) > 0, KsiR) > and all R> 0. 

Note that (/12) and (/13) are satisfied, e.g., by 6(a) = A(a) = e"/"^, a > 0, with r > 0, which is one 
of the choices of A and 6 in [31 [31 [HI [IHl [IS]- Also note that the function D in (^,5) may be such 
that D{0) = 0, that is, we may allow for a degeneracy of the diffusion coefficient D{r) at r = 0. 
For example, D{r) := D^r^ with > and 6* > 1 (or = 0) satisfies {h^). Following [i], the 
function E{r) := D'{r) = ODqt^^^ fulfills assumption (/ig) in this case. 

Definition 2.1. Suppose (hi) — (Hq). A (global) weak solution to (|l.ip - (|1.6p is a pair of non- 
negative functions (u, v) possessing, for each T > 0, the regularity 

u G Loo(0, T; Li((0, 00) x n; 6(a)dada;)) , v G C^[0, T];L^(n)) , 



A := / A(a)u(.,a,.)daGLoo((0,r)xr!), ^ L2iO,T;W^{n)) , j-1,2, 

Jo 

and satisfying v{0) — v'^ , 

/•oo 

dtv — {g{v) — S,{v)) V + / 6(a) /i(a) u(., a, .) da a.e. m (0,T)xO, 



and 



{dt(p + da^p — tJ. (p) u dadxdt + / / if{t,0,x) £,{v{t,x)) v{t,x) dxdt 



Jo. Jo Jo Jn 

00 



/> /'OO />! /> /'OO 

+ / / (p{0,a,x) u^{a,x) dadx + / / / '^xf D{A) u dadxdi 
Jn Jo Jo Jn Jo 

-Li { w ^'f''^' - fw ^"^-'^'l ^ 

/or any test function ip G C^([0,T) x (0, 00) x £7) with compact support and d,^(p{t,a,x) — for 
it,a,x) G (0,T) X (0,00) X 
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Observe that the regularity required on u, v, and A along with assumptions (/13), (^5), and 
(/ig) ensure that all the terms in the weak formulation for (jl.ip are meaningful. In particular, 
E{K,v)/C,2(K) and £''(A)/C((A) are both bounded by (/15) and {h^) while (/13) implies that 



u Vxf da 



AeLo,((0,T)xO), 



Let us emphasize here once more that only the averages of u with respect to age have the needed 
intcgrability properties for the weak formulation to make sense. A related situation is encountered 
in the theory of renormalized solutions for the Boltzmann equation and parabolic equations 
(see, e.g., [HISIIS] and the references therein). 

Our main result is then the following: 

Theorem 2.2. Suppose {hi) — (hg) and let p > N. Then, given any non-negative initial values 
{u'^,v^) satisfying 

w° e Li((0,oo) X n;b{a)dadx) H Li{{0,oo),W^{n); X{a)da) , t>° e W^i^) , (2.1) 
w° lnM° e ii((0, 00) X n; A(a)dadx) , (2.2) 
there exists a global weak solution to (ll.l[) - (jl.6p . 

The regularity assumptions on the data D, b, A, and the initial data u", v'^ could be weakened, 
see Remark 14.71 

3. A Regularized Problem 

The basic idea to handle age structure is to discretize equation p.ip with respect to the age 
variable a G (0, 00). To that end we fix / G N and a G (0, 1). Let then 

{bi)i<i<i+i , {K)i<i<i+i , and (/^i)i<,;</+i 

be non-negative numbers such that 

k>l, A, > ^ > , 

< b* < Bb, , A* < L A, , ^J.^ <M , |Ji^h < I3K < h , 

for some positive numbers ^, L, M, and /3, where we put 

,» bi+i — bi Xi+i - Xi . 

:= , A,; := , 2 = 1,...,/. 

a a 

Moreover, let O £ C°°(M) be a cut-off function satisfying 

e>0, e' < , e(r) = 1 for r<l/2 and e(r) = for r>l. (3.2) 
Finally, let ^, g, D, and E be functions satisfying (hi), (/15), (h^) together with 

D{r)>do>0, {l + s)£,{s) + E{r,s) <E, (r, s) G M^ (3.3) 
for some constants do > and S > 0. 

We then look for a solution {ui, . . . , uj, A, v) to the approximating problem 

dtUi + —{ui - Ui_i) = div:r (-D(A)Vj;Ui + UiQ{a'^Ui) E{A,v) V^A)) - fi^Ui , (3.4) 



(3.1) 



aA = di 



L'(A) + ^ aAiUi£:(A, w) 



VxA + Xi£^{v)v + ^ q;(A* - HiXi)Ui - A/+1 m , 

J i=l 

(3.5) 

dtv = alS.xV -t- [g{v) - £,{v))v + a ^ bifiiUi (3.6) 



I 
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for i = 1, . . . , / and {t, x) S (0, oo) x fi, where 

uo ■■= £.{v)v , (3.7) 

and subject to 

d^Ui = d^A :^d^v^O , (t, x) E (0, oo)xdn , (3.8) 

and 

u,{0,x) = u°{x) , A{0,x) = A°{x) -.^ a^Xju"j{x) , v{0,x) = v"{x) , xeil, (3.9) 



for i — 1, . . . , /. We will first show that 
some uniform bounds on this solution. 



possesses a classical solution and then derive 



3.1. Global Existence. In this subsection we prove the global well-posedness of p.4 p- p.9p . More 
precisely, we have the following result. 

Proposition 3.1. Let p > N , a <E {0,1), and consider {uf, ... ,u'j,v°) e Wp{n,R^^^) with 

< u1(x) < , < v"(x) , xen . (3.10) 
4 

Then there exists a unique classical solution 

{ui,...,ui,A,v) e C([0,oo) X n,R^+2) nCi'2((0,oo) x Cl,R^+^) 

to problem p.4p - p.9p such that Ui{t, x) > and v{t, a;) > /or i = 1, . . . , / and (t, x) G [0, oo) x fi. 
Moreover, setting 



t* :=sup<i>0; max sup ||Mi(r)||oo < i >0 



l<i<I 



-G[0,i] 



2a2 



have 



A{t,x) = a^\iUi{t,x) , {t,x) e [0,t*) X n 



(3.11) 



(3.12) 



Proof. We fix rj such that 



< 77 < 



do 



a ||i?||oo / max Xi 
i<i<i 



and put Do := (-?/, 00)^+2. Moreover, we define a := (a"'") e C'^ {Do, C{M.'+^)) by 



/ Diyi+,) 

Diyi+,) 



a(y) 



yiQ{a^yi)E{yi+i,yi^ 



\ 



V 







D{yi+i) yiQ{a'^yi)E{yi+i,yi+2) ■■ 

I 

D{yi+i) + aJ2 ^iE{yi+i,yi+2)yi 

4=1 

a J 



for y ^ (yi, . . . , yi+2) e -Do- We next set Oj^kiy) := a{y)5j,k for I < j,k < N and y G Dq, and we 
introduce the operators 



N 



N 



A{y)z -.^ - ^ dj{aj^k{y)dkz) , B{y)z -.^ ^ a.j^k{y)dkz 

1 j,/c — 1 
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for z — (zi, . . . , 27+2) and the function 



<?ri</+2 



-A^ij/i — (yi - ^{yi+2)yi+2) 

a 



-M22/2 — {v2-yi) 

a 



1 



AiC(2//+2)y/+2 + a^(A* - ^.i\i)y, - Xi+iyi 



V 



(5 - Oiyi+2)yi+2 + a^h 



'i^J■iyi 



With these notations, an abstract formulation of (|3.4p - (|3.9p reads 

dtZ + Aiz)z = fiz) , 
B{z)z = , 
z(0) = K,...,u?,AO,t;0) . 



Clearly, owing to (|3.3p and the choice of 77, the eigenvalues of a{y) are positive for each y e Dq, and 
the boundary- value operator {A,B) is of separated divergence form in the sense of [fl Ex. 4.3(e)]. 
Consequently, the boundary-value operator {A, B) is normally elliptic. It then follows from [H 
Thm.14.4, Thm.14.6] that (|3.4p - p.9p has a unique maximal classical solution 



z = {ui,...,ui,A,v) e C([0,t+) X n,Do) nCi'2((0,i+) X n,R'+^) , 

where i+ G (0, cxo] denotes the maximal time of existence. Observe that a^'"(0,y2, 
forn e {2, . . . , / -f 2} and /^(O, y2, ■ ■ • , yi+2) > by (/ii). Therefore, [H Thm.15.1] ensures that 
ui{t,x) > for {t,x) G [0, i+) X Q. Now, U2 solves 



,yi+2) 







atM2 - div^ {D{A)\/^U2 + U2e{a^U2)E{A, v)\/^A) 



M2 



U2 = —Ui > 

a 



with a non-negative initial condition, which readily entails that U2{t,x) > for {t,x) G [0,t~'') x J7 
by the comparison principle. Proceeding by induction, we obtain in a similar way that Ui{t, x) > 
for {t, x) G [0, i+) xH and i G {1, . . . , /}. The same argument gives v{t, x) > for (t, a;) G [0, t^)xO,. 
We next show that — 00. To that end we define the parabolic operator Ci by 

Ciw dtw - div^{D{A)Vx'w) - (a^'Uie'(a^wi) + e(a^wi))£'(A, 'y)Va;A • V^rW 



div^ (£'(A, v)V^A) w Q{a'^w) 



Setting k{t) + Et/a > Yjo? for t > 0, we infer from (/ii), (|3.3p . and the properties of 

that C.\k > k' —^{v)v/a > 0. Therefore, ui{t,x) < k{t) for {t,x) G [0, i+) x fi by the comparison 
principle since we assumed that ui(0, x) < IjoP' for a; G 17. Furthermore, we have C2U2 — 0, where 
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£2 is the parabolic operator defined by 

Liw := dtw - diV:r(-D(A)V:rw) - {a^ {o? u-^) + Q{o?-U2))E{k,v)^ ^k-\I 



Ml 

a 



Owing to the previous bound on ui and the properties of 0, we have £2^ > k' + (k — ui)/a > 0. 
Using again the comparison principle we conclude that U2{t,x) < k{t) for (t, x) G [0, t^) x Cl. 
Proceeding analogously for Uj, i G {3, . . . , /}, we derive that 

0<Ui{t,x) <k{t) , (i,a;) e [0,i+) X , {!,...,/}. (3.13) 

Now, by [hi), p.6p . and (|3.13p we clearly have 

I 

dtv - alS^xV < II5II00W + ak{t)'^bi fit 

i=l 

from which we deduce that 

/ /■* 

< ^(i,^) < ||i;°||ooe"»"°^* + a Vfo.Aii / A:(s) eUsH-^'^^'Ms , (i, x) G [0, t+) x 17 . (3.14) 
Finally, by (/ii) and (|3.13p we have 

-\i+ik{t) - q;^(|A*| + fii\i)k{t) < Xi£,{v)v + ^a(A* - ^iiXi)u., - A/+i u/ 

i=l i=l 

I 

< Ai5 + aY,\X*\k{t), 



i=l 



so the comparison principle applied to (j3.5p warrants that 



m,x)\ < ||A"|U + AiSi + 



A/+1 + (|A*| + fiiXi 



1=1 



fc(s) ds (3.15) 



for {t,x) G [0,f+) x ft. Thanks to (|3.13p . (|3.14p . p.lSp . and the upper triangular structure of the 
diffusion matrix a, we are in a position to apply [2 Thm.15.5] and conclude that indeed t'^ = 00. 

It then remains to check (I3.12p . First note that > due to (|3.10p and the continuity of 
[ui, . . . ,ui). Next, setting P := aJ2i=i \ui, it follows from (|3.4p that P solves 

dtP + ^Xi{ui - u,^i) = div^r {D{h.)\/^P + a^XiU,Q{.a^Ui)E{A,v)\I^K \ -a^X^^,u, . 

j=l \ 2=1 / i=l 

On one hand, we clearly have 

/ / /-I / 

Xi{u^ - Ui^i) = ^ XiUi - ^ Aj+iUj = -XiS,{v)v - a ^ X*Ui + A/+iU/ 



2=1 i—1 i—0 i—1 



by p.7p and the definition of A*. On the other hand, if i G [0,t*), then Q{a'^Ut{t, x)) = 1 for x G 
and so C(J2i=i XiUiQ{a'^Ui) = P in [0, <*) x Cl. Consequently, P solves the same initial-boundary 
value problem (|3.5p . (|3.8p . (13. 9p as A in (0,t*) x il. The uniqueness of classical solutions to this 
problem guarantees that p.l2p holds true. Thus the proof of Proposition [57T] is complete. □ 
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3.2. Uniform Bounds. The aim of this subsection is to derive some uniform bounds on the solu- 
tion obtained in the previous subsection. Given p > N and (u?, . . . , Uj, w°) G Wp (O, R^+^) obeying 
p.lOp . let (ui, . . . , ui, A, v) denote the classical solution to p.4p - (|3.9p provided by Proposition l3.1l 
We fix a constant Kq such that 

[ faV6^w°(x)+w*'(a;)+aVA, [u^(a;)(lnu^(a;)-l) + l] ) dx 

\ ^=l ^=l / 



+6i + Ai + ||A"||^ + ||t;^||^ < K„. (3.16) 

In the following, c and Cj, j > 1, are generic constants that may differ from place to place and 
depend on £, B, L, M, /3 in p.ip . ||g||oo, and Kq, but not on /, D, E, a £ (0,1), do, and S in 
p.3[) . Dependence on additional variables will be indicated explicitly. 



We start with an ii-estimate: 
Lemma 3.2. For T > 0, we have 



J (^aJ2b^u^{t,x) +v(t,x)^ dx < ci(T) , t e [0,T] . (3.17) 

Proof. Multiplying p.4p by abi, summing with respect to i, and integrating with respect to a;, we 
obtain from (|3.7p and p.Sp the identity 

— / a biUi dx — 6i / ^{v)v dx — f a b* utdx ^ — f (a bi fii Ui + u/ | dx . 
Jo .^^ Jo Jo -^^ Jj2 \ -^^^ y 

Integrating ()3.6|) with respect to x and adding the result to the above identity gives 

— / { a'S^'biUi + v\ dx < [ (3(7;) — (1 — 6i)^(u))udx + / a'S~\b*Uidx 
dt Jn \ J Jn Jn 









/ w dx + 


i sup 




In 


Vi<j</ 



due to {hi), whence the claim from (13. ip . □ 

We next improve the previous Li-estimate. An appropriate choice of the sequence (ryi)i>i con- 
sidered in the forthcoming lemma will allow us to control the tail of the approximating sequence 
(ui) in Li((0,r) X (0,00) X r2; 6(a)dtdadx) later on. Indeed, such a property in turn will guar- 
antee one of the two conditions required for the Li-weak compactness of the approximation (see 
Section 14.31 below) . 

Lemma 3.3. Let (rji) £ [0, 1]^+^ be such that rji — and rji < rji+i for i = 1, . . . ,/. Then, for 
T>0, 

/ aVr?,6,M,(i,x)dx < e^* /" a V 7y, 6, u^(x) dx jry* L C2(T) , te[0,T], 
where rj* := (77.^+1 - r]i)/a for i = 1, . . . ,1 and \r]*\^ := maxi<i</ \r]*\. 

Proof. Multiplying p.4p by otrjibi, summing with respect to i, and integrating with respect to x, 
we obtain from (|3.7p and (j3.8p the inequality 



37 / a -qihui dx - [ aS^ (bi+i 
diJo ^ Jn tt\ 



f Oj Ui dx < . 



a 
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Owing to (|3.1|) we have — bi + a b* < {1 + aB) bi, so that 

^ / a riihui dx < / « ((1 + "-B) r]* + ij^ B) bi Ui dx 
di Jn ,_i Jn 



i=i i=i 



< B / aVry,&,u,dx+(l + S) |ry*|^ ci(T) 



i=l 

the last inequahty being a consequence of p.l7p . The claim then follows by integration. □ 
We next derive Loo-estimates on A and v. To this end recall that > was defined in (|3.f f p . 
Lemma 3.4. For T > 0, we have 

||A(<)|U + lk(t)||oo <c3(r) , i e [o,r]n[o,C) • (3.18) 

Proof. We first infer from (hi), ^Mj, and ((31^ that 

dtv~aA^v<\\g\\^v + P A in (0,C)xf7 
with homogeneous Neumann boundary conditions. The comparison principle readily entails that 

\\vmoo<(\\v"\\^+P l^\\Ais)\\^ds^ eNI-*, te[0,C). 

Next observe that, by (hi), (|3.ip . (|3.5p . and (|3.12p . we have 

dtA~diY,i[D{A)+AE{A,v)]V^A)<\i\\g\\^v + LA in iO,t*Jxn. 
Using the comparison principle once more, we deduce that 

||A(t)||oo< f||A°|L + Ai||.9|U f \\vis)\\^ ds) tG[0,C)- 



Consequently, 



l|A(i)||oo + lk(t)||oo<c(T) (^1 + i\\A{s)\\^ + \\v{.s)\Udsj , ie[o,r]n[o,t:), 

from which the claim then follows. □ 

As a consequence of the preceding lemma we obtain a lower bound for i* . 
Corollary 3.5. Consider T > 0. If a < £/{4c3{T)), then t*^ > T. 

Proof. We consider T > and a < £/(4c3(T)). Assume for contradiction that < T. Then it 
follows from p.ip . (|3.12p . p.lSp . the non-negativity and continuity of Ui, and the choice of a that 

a£u,{tl,x) <A{t*^,x) <C3{T) <^ , xeQ, t = !,...,!, 

whence Ui{t*^,x) < l/(4a'^) for a; e il and i = Together with the continuity of the Ui's, 

this contradicts the definition (|3.1ip of t* . □ 

Next we establish some bounds on (ui) that will guarantee its local weak compactness in Li 
and the strong compactness with respect to space and time of its averages with respect to the age 
variable. As already mentioned, this approach is inspired by the existence proof of renormalized 
solutions to the Boltzmann equation [7 which makes use of velocity averaging results, see, e.g., 
[T2] . [TTl Chapter 5], and the references therein. 
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Lemma 3.6. For any T > 0, 

y^aXi I <t>{ui{t,x))<lx < Ci{T), (3.19) 
Jn 

aJ2X^D{A)\V^uy^\^ + E{A,v)\V^A\^\ dxds < Ci{T) , (3.20) 



Jn 



for t e [0, T] n [0, i* ), where (j}{r) r(lnr - 1) + 1 /or r > and 0(0) := 1. In addition, 

(|V,Ci(A)|' + |V,C2(A)|') AxAs < C5{T) , te[0,T]n [0,0 , (3.21) 



/o Jn 

where Ci and C,2 are defined in (h^) and {he), respectively. 

Proof. Multiplying p.4p by aXilnUi, summing the resulting equations with respect to i, and 
integrating over f2, we obtain 

/ I ^ I 



^ J Xi 4>{ui) dx - ^ y Ai (wi - Ui-i) Inw^ da; - "^XJ j ^'^^^ 



V u |2 

" + £;(A, v) \7^u^ ■ A !> dx 



Taking into account that 

— rlnr < (/)(r) + r , r >0 , 

and that 

{ui - Uj_i) Inui > (j){u,) - (j){ui^i) 
due to the convexity of (f>, we derive from (|3.7p and ()3.12p 



^ ^ a Aj y 0(Mi) < - ^ Ai y (0(ui) - ^(ui-i)) dx + ^ y a Ai /ij (0(wi) + Wi) dx 

- / |^aA,i?(A)J 
U=i 



1^-^^ + E{A,v) |V,A|' dx 



< Ai / (/)(^(w)u) dx — A/+1 / (/)(M/)da; + aTJ / X*(j){ui)dx 
Jn Jn Jn 

+ aAj^j ((/)(mj) dx 

-/ |^aA,i?(A)^^^ + i;(A,i;) |V,A|'| dx . 

Recalling that (p > 0, {hi), (|3.ip . and (|3.18p we deduce that 

VaAi / 0(uOdx < c(T) + a (L + A/) V / Ai0(w,)dx 
d* i=i ■''fJ .=1 

-y |^aA.Z?(A)r^2^ + £;(A,z;) |V,A|'| dx , 
from which ([XTO)) and (P?^ follow. 
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Next, by the Cauchy-Schwarz inequality and (j3.12p we have 



|V,A|<2 (j2aK\y,uy^\'^ 



1/2 



From this, (/15), and (|3.20p . we deduce that 

/ / |V,Ci(A)|' dxds = / / 
Jo JQ Jo Jn 

< 4 

while (he), |3J8| . and ((3:20)) imply that 



|V^C2(A)| dxds < 

/o Jn 

thus completing the proof. 



^(A) 
A 



|V2;A| dxds 
I 



Jn 



D{A) ^aA^IV^u^y dxds < C4(r) , 



i=l 



E{A,v) 



Jn i^2{c3{T)) 



«2(c3(r)) ' 



□ 



Lemma 3.7. Consider (xi) £ ^^'^^ with xi+i = '^'^'^ P'ut :~ ctX]i=iXi^i '^^^ %* := 
iXi+i ~ Xi)/'^ for i — 1, . . . , /. Then, for t G [0, T] n [0, i* ) T > 0, we have 



{\\\/^V^{M^)\\l + \\dtV^{M^)\\^wk^^^n)y) < cg (T, |x|oo, |x*|oo, G{T)) 
where 

G{T) := ||£)||i^(o,c3(T)) + II-S|Il^((o.c3(t))2) , 
anii g C^(R) is anj/ function satisfying T>^{Q) ~ ^?^(0) — and 



(3.22) 



< T>''{r) < D 



Ixic 



r e . 



Proof. First observe that the Cauchy-Schwarz inequality together with (|3.ip and (j3.12|) give that 



1/2 



1/2 



2 a|xU > 1/2 



1/2 



< 



2 Ixic 



1/2 



Ai/2 ^aA. 



1/2 



Consequently, by (|3TT8| and 



f I D{K) |V,A^J' dxds < twAmoo f i^aKDiA) 

Jo Jn Jo Jn j^]^ 



< c{T) \xL 



1/2 



Owing to (|3.ip we have 



lA^xl < Ixloo A/£ 



and the monotonicity of D warrants that 



D 



Ixic 



dr < \Mx\ D 



Ixloo 



dxds 
(3.23) 

(3.24) 

<\M^\D{K)<''^kD{K). (3.25) 
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We then infer from (jXTBl) . ^M), and (jX^ that 



Jn 



£2 



D{Af \V^M^\^ dxds 



Jn 

2 



< <T)\x\L, 



whence the first part of (|3.22p . 



Next, we infer from (|3.4p that Aix solves 



dtMx = div^{D{A)\7^Mx + MxE{A,v)\/^A)+Mx*-xf,+Xi^{v)v in (0,i*)xfl 

with homogeneous Neumann boundary conditions. Consider ip £ W]ij^i{fl). Multiplying the above 
equation by ip 27^(A^^) and integrating over Q give 



where 



Fi := 
F2 := 



/o Jn 



{D{A) \7^Mx + E{A, v) V^A) V'^iM^) V^ip dxds , 

f f {D(A) V^M^ + E{A, v) V,A) P^(Xx) V.A^x V dxds , 
Jo Jn 

(7Wx*-xA' + XI I') %{Mx) V dxds . 

that 

1/2 



(3.26) 



"'ii 



It first follows from (jXTS)) . (IX^ . (jX^ . and that 



\Fx\ < ^ ||i?(A)i/2 V^A^^l^ (^^^(^) ^x(-^x)' |Vx<Ppdx) ds 



+ ^ ||f(A,z;)1/2 v,a||_^ ^ F(A,t.) Ml V'^iM^y \S7^p\^dx^ ds 

< C{T, \x\oo) ||V:r</?||A'+l • 

Next, p.24p . the assumption on V^, and the monotonicity of D yield 



\ Ixloo 



<D[A). 



We then similarly infer from p.lSp . p.20p . (|3.23p . and the continuous embedding of Wl^j_-^{Vl) in 
Loo(f^) that 



\F2\ < \W\\oc I I D{A) V'^iM^) iV^M^f dxds 
Jo Jn 



E{A,vf'^ V,a||_^ (^j E{A,v) Ml V'^^iM^f \V^Mx\^dx^ ' ds 



< cM^.^^^^n) / DiAy \V,Mx\' dxds 



Jn 



+ cM\w^^^^in) -^x ^(^)' l^xM^l^dxds^ 

< c{T, Ixloo) \Mw^^ (n) ■ 



1/2 
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Finally, we deduce from (hi), ([Sj]), ([XTS]) . and ([X^ that 

< c(r, Ixloo, Ixloo) ||<^||w^^^(n) ■ 

Combining the above three estimates for Fi, F2, and F3 with p.26p leads to the assertion by a 
duality argument. □ 

We finish off this section with an estimate on AxV. 

Lemma 3.8. For any T > 0, we have 

f WA^llds < a\\V,v"\\l+c{T) , i e [0, T] n [0, C) • (3.27) 
Jo 

Proof. We multiply p.6p by —aA^v, integrate over Q, and use {hi), (|3.ip . and (|3.18p to obtain 

f ^l|V.z;||^ + a2||A,«||2 < allfflU IIHI2 ||A.«||2 + a/3||A||2 ||A,^;||2 

2 

< ^||A,Hl2 + c(T) 

in [0, T] n [0, <* ). Integrating with respect to time yields (P??7)) . □ 

4. Passing to the Limit 

We are now in a position to prove Theorem l2.2l We thus construct a weak solution to (|l.ip - p.6p 
in the sense of Definition 12.11 by using a compactness argument for the solution to the regularized 
problem (|3.4p - (|3.9p . First we demonstrate how we set up equations 



4.1. Approximation. Suppose hypotheses (/ii) — (/ig)- Choose a G (0,1) and put /(a) := [l/a^]. 
We then set 

Da{r) -.^ D{r) +a , reM, 

and note that satisfies p.3p with dp — a. By classical approximation arguments, we also 
construct a non-negative function G C'^(M^) fl Loo(R^) such that Ea{r, s) = E{r, s) for (r, s) G 
[0, 1/a]^ and a non-negative function G C^(R) fl Loo(R) such that s (1 + s)(^a{s) belongs to 
ioo(K) and ^q(s) = ^(s) for s £ [0, l/a]. Furthermore, we set, for i = 1,. . . ,I{a) + 1, 

K,a ■= — X{a)da , bi^a ■= — 6(a) da , fii^a := - / n{a) da 



and, for i = 1, . . . , /(a). 

Observe then that (^12) — (/14) imply the validity of (|3.ip with _B = Bq, = -^0, /? = /3o(l + -Sq), 
and ^ = ^0- Indeed, hypothesis (ft-2) ensures that 

0<b:^ = ir ^(" + ")-^(") da<i?ok., . = l,...,/(a). 

One shows A*^ < LgXi ^ analogously using (^,3) which also gives A^^^ > Iq. Finally observe that 
(hi) warrants 

fJ-i,abi,a = ^ / M(a) / 6(2;)dzda < -| / A(a) / — - dzda 

<^ J(i-1)Q J(i-1)Q " J{i-l)a J{i-l)a Oyo.) 
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and therefore Hi^a bi^a < Po{Bq + 1) Xi^a owing to 



b(z) b(ia) „ , / ^ 

5(a) b{{i - l)a)) 



since b is non-decreasing. Given p > N and any non- negative initial values {u^,v'^) satisfying (j2.1l 
and (|2.2p . let (ui^q, . . . , u/(a),Q, Aq,, denote the classical solution to ()3.4|) - (j3.9p with / — I{a) 
{Xt,bi,fi,) = {Xi,a,bt^a, ^J■i,a), and where {D,E,^) are replaced by {Da,Ea,S,a) and u° by 



a(^) ■= ~ / u°(a, x) da , X G , z = 1, . . . , /(a) . 



(i-l)Q 



Note that we may assume without loss of generality that ||u^^||oo < l/4a^ by making a smaller if 
necessary. We first collect some properties of (u^^^, A°). 



Lemma 4.1. For a > small enough, we have 

n /(«) 

/ a^(&,,„ <„(x)+A,,„ </)«„(x))) dx+llA^II^ <C7, 

where (j){r) := r (Inr — 1) + 1 for r > and (t>{0) :— 1. 
Proof. By (/la), 

I- Ha) ^ „ I{a) „ia i-ia 

/ ay^bi^a a(^) — — 2, / M°(a,x) / b{z) dzdadx 

in ' ' ^ -^^i^ J{i-l)a J{i-l)a 

< i y b{a)u^{a,x) / ^i^±^dzdadx 

a Jn ~[ J{i-i)a J{i-i)a b{a) 



< (l + Bo) / / 6(a)u"(a,a;) dada;. 
Jn Jo 
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Next, (h^), (|2.2p . and Jensen's inequality entail that 



4=1 



= I Y I I Hz) dz] (pi- f u"{a,x) da] dx 
^ - / y / Hz) / (f>{u°{a,x)) dadzdx 

^ Jn~l J(i~l)a J{i-l)a 

= — / 2, / / A(z) dzdadx 

^ J ~[ J (i-l)a J{i-l)a 

+ - y2 <P{u°{a,x)) / A(z) dzdada; 

1 r ^'-"^ ( r \ 

^- ill A(a)0(M°(a,a;)) / {1 + Lo{a ~ z)) dz ] dadx 

^ Jn'~[ J(i-l)a \J{i-l)a J 

+ - j Ylj^ A(a)0(u"(a,a:)) iyj {1 + Lo{z - a)) dzj dadx 

I Ha) „ia 

<(l + io) y2 X{a) (l){u"{a,x)) dadx 

Jn . , J(i-l)a 



<(l + io) / / A(a) ^(^"(a,^)) dado: 
Jn Jo 



Finally, as above we deduce 



0<A°(x) = - y / u"{a,x) / A(z)dz+ / A(z) dz da 

" ~! J(i-l)a \J(i-l)a J a / 



< (l + Lo) / A (a) u" (a, a;) da 



and the right-hand side of the above inequality belongs to Loo{i^) as a consequence of ()2.ip and 
the continuous embedding of Wp{il,) in Loo(f^) (recall that p > N). □ 

Now we introduce 

/(a) /(a) 

\a{a) := y Ai,al((i_i)Q,^iQ](a) , 6a(a) := ^ &i,Ql((i-i)Q,m] (a) , 



i=l 4=1 
/(a) 



A*Q(a) := y A^i,al((i-i)a,iQ](a) , 

4=1 

and 



(i, a, x) := y W4,a(^, a^) (a) 



4=1 
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for {t,a,x) e [0,00) X [0,00) X Q. Let 

:= sup < t > ; max sup ||ui.Q(r)||oo < 77^ i > , 
[ i<'<-f re[o,t] 2a ) 

be defined as in (|3.1ip . We first estabfish that, as expected, is a weak solution to an approxi- 
mation of the original problem. 

Lemma 4.2. Ifcp^C^ ([0, 00) x f2) is such that supp (p C [0,R]xCl for some R> and al(a) > R, 
then 



d f°° f Iff 

— ipuadadx — — / (p£_a{va)vadxda 

Jo Jn a Jo 

' ip{a + a) — ip{a) 



Jn 



— (^(a)/iQ,(a) I Ua{a) dxda 



Aq) • Vx'pdxda 

Jo Jn 



forteiO,tl). 

Proof. Let ip £ C^([0, 00) x f]) with suppiy3 C [0, R] x and put 

Vi,a{x) := — / ip{a, x)da , z = 1, . . . , /(a) + 1 , x £ . 

" J(i-l)a 

We infer from ((331) and jlJl) that, for t e (0, i;^), 

/(a) 



17 / / Mq dadx = ^ / a Ui,Q da; 
Jo dt ^ 



i=l 
/(a) 



dx 



a 

i=l -^Sl 

= / ¥'l.QCQ(Wa)Wada; + / {(pi+l^a - '^i,a)ui^adx - / 93/(q) + 1,q a dx 

Jn Jn 

— / / 93 /ifi Uq dxda — / / (£'Q,(Act)V Aq) • Vx</3dxda . 

Jo Jn Jo Jn 

Noticing that al{a) > R implies ^i(a)+i,a = 0, the assertion follows. □ 

4.2. Compactness estimates. Our aim is then to pass to the limit as a in the identity 
stated in the previous lemma. We thus need to provide some compactness for (ua), (Aq), and 
(va), a first step being the derivation of suitable estimates. 

Lemma 4.3. Let T > 0. Then, for a small enough (depending on T), A > 4, and t e [0,T], the 
following estimates are valid: 



Jn 



\b{a) Ua{t, a, x) + A(a) (f>(^Ua{t, a, x))] dxda < ci(T) + C4(T) , (4-1) 
IIAqWII WI|oo<C3(T), (4.2) 

lo In + |V.C2(Aq)|' + a |V,Aq|') dxds < cs{T) , (4.3) 
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/ / b{a)uait,a,x)dadx < csiT) iff b{a)u°{a,x)dadx + ^] . (4.4) 
Jo J A \Jn J A/ 4 ^ J 

In addition, for any x G ^""^([0, oo)) with compact support, the sequence (M^^a) defined by 

POO 

Mx,a{t,x) := I x{a) Ua{t,a,x) da , {t, x) £ {0,t*^) x fl , 
Jo 

is such that 

ll-^x,a(OIL+ / {\\^MMx,o.)\\l + \\dtV^iM^..c.)\\^wk^,m') ds<c9(T,||x|ki(o,oo)) (4.5) 
for t e [0,T] n [0,tl), where the function e C'^{M.) is defined by V^{0) = P^(0) := and 

Proof. Owing to Lemma HJI and assumption (ft.4), the condition (|3.16p is fulfilled with 

Ko C7 + (l + |l}|)||w°||oo + (l+/3o)K2) , 

which clearly does not depend on the approximation parameter a. 

Let T > 0. According to Corollarv l3.5[ we may choose a small enough (depending on T) such 
that t* > T. We can also assume that a satisfies a < l/cs{T), the constant C3(T) stemming from 
in Lemma [331 

Observe first that (|4.H) and ()4.2|) are immediate consequences of Lemmata 13.21 13.41 and 13.61 A 
useful consequence of (|4.2p and the choice a < l/c3{T) is that ^a{va) — ^{va) and Ea{Aa,Va) = 
E{Ka,Va)- We then infer from ([^^ that 

^ ^ (|V,Ci,a(A„)f + |V,C2(A„)|') dxds < C5(T) , te[0,T], 

with Ciair) ■■= ((Z?(r) +a)/r)i/2 and Ci,a(0) 0. Clearly, C'lA^c) > C((Aa) and also 

ai/2 ai/2 



by (|4.2p . Collecting the above information allows us to conclude that (|4.3p holds true. Also, 
introducing 

Xi-=- / x(a) da and x* ■= , « > 1 , 

a J(i-l)a a 

for X G C^([0, 00)) with compact support, we readily see that |x|oo + |x*|oo < 11x11 wj^(o, 00) and 
(|4.5p follows at once from Lemma [5T71 (|4.2p . and the inequality Mx,a < ||x||ooAq/^. 

Finally, let rj G C°°(K.) be a fixed non-decreasing function such that 77(a) = for a < 1/2 and 
r](a) = 1 for a > 1. For A > i and « G N \ {0}, wc put rji := rj{ia/A). Then 771 = and {rii)i>i 
is clearly a non-decreasing sequence with < rj* := (ryi+i — rji)/a < ||ry'||oo/^ for i > 1. We then 
infer from Lemma [5T51 that . for t G [0, T], 



/(a) „ /(q) 



/ " '?i ^'i,a(*: a^) da; < e^"* / ay2mbi,au'^ ^{x) dx 
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The properties of ba and rj imply that 



b{a) Ua{t, a, x) dadx < / / ba{a) Ua{t, a, x) dadx 
in J A JnJA 



< 



< e^°* / / 5„(a)M°(a,x)dada: + C2(r) , 



lnJ{A-2)/2 ^ 

whence (|4.4p by (/12), the latter guaranteeing that ba{a) < (1 + Bo)b{a) for a > 0. □ 



Remark 4.4. Owing to the superlinearity of at infinity, the estimate (|4.1[) warrants the weak 
compactness of (ua) in Li((0,T) x (0,^) x fl) for T > and A > and strongly relies on the 
assumed positivity of X. Therefore, if X would vanish on some interval (0,ao) with Oq > (as, e.g., 
in [21 131 [HI HI] ^, the restriction of to the set (0, T) x (0, ap) x fl is only weakly-* compact in 
the space of bounded measures. In that case the passage to the limit performed in the next section 
might be more delicate. 

4.3. Proof of Theorem [2721 

Proof. We fix T > 0. Due to the positive lower bound on Aq, (|4.ip . and (14. 4p . we may apply the 
Dunford-Pettis theorem [51 IV. 8] to conclude that there are a sequence {ak)k>i with offc — > and 
a non-negative function u G Li((0, T) x (0, 00) x il; b{a)dtdadx) such that 

Ua^ u in Li((0, T) x (0, 00) x Q.; b{a)dtdadx) . (4.6) 

We may also assume that, for each fc > 1, a*: is small enough such that i*^ > T and c^{T) < I/at. 
Next, given any x ^ C^([0, 00)) with compact support, we readily deduce from (14. 6p and the 
positivity and unboundedness of b that 

M^.c.,~-M^ in ii((0,T)xn), (4.7) 

where A^^^q,^ is defined in Lemma l4.3l and 



Mx(.t,x):^ / x{o) u{t,a,x) da , (t, x) e (0, T) x . 
Jo 

Furthermore, owing to (14. 5|) . we may apply ,19, Corollary 4] to conclude that (2?^ (-^x,"fc))fc 
relatively (strongly) compact in L2((0, T) x $7), hence converges also a.e. (after a possible extraction 
of a further subsequence). This property, the strict monotonicity of T)^, and (14. 5p then imply that 

A^x,"fc — ' in ig((0,r)xl7) for any qe[l,oo). (4.8) 

We next claim that 

— * A:=7Wa in Lq{{0,T)x^) and a.e. in {0,T)xVl, (4.9) 

Xbc.P^^.a. ^ -^ftM in ig((0,r) X r!), (4.10) 

for any q G [1, 00). Indeed, let •§ £ C°°(K) be a smooth and non-increasing cut-off function satisfying 
d{a) = 1 if a < 1/2 and i9(a) = if a > 1. For A > 1 and a > 0, we put ^^^(a) 19(0/2^). We 
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infer from (/14) and (|4.4p that 



r I |A«, - A| Axdt = f I I 
Jq Jn Jo Jn Jo 

T 



< 



Jn 

Po 

T 




T 

Jn J A 



< 



Jn 



+ c{T) 



A(a) (uq^ — u) da dxdt 
A(a) i?yi(a) (wqj, — u) da da;di 
6(a) (uqj. + u) dadxdt 
\M\i}A,ak - dxdt 

b{a) (w° + u) dadxdt + 



JnJA/4 



A 



whence 



Urn sup [ I - K\ dxdt < c{T) ([[ [ ' 
fe^oo Jo Jn \Jo Jn J A/4 



6(a) (m° + u) dadxdt + 



A 



by (|4.8p . Letting A ^ 00 completes the proof of (|4.9p for q = 1. The extension to g G (1, 00) next 
follows from (14. 2p by interpolation. The proof of (|4.10p is similar and uses additionally (/i2)- 

A further consequence of (|4.8p is that, for any f e C^([0, T) x [0, 00) x fj) with compact support 
and (7 G [1, 00), we have 



(4.11) 



(^(t, a, x) (i, a, x) da and A4ip{t,x) := / (^(i, a, x) u(t, a, x) da 



with the notations 



for (ijx) G (0,T) X i7. Indeed, we argue as in [71 Section IV] and first note that (|4.1ip readily 
follows from (j4.8p if there are an integer J > 1 and functions ('0j)i<j<J in C^([0,T) x fj) and 
{Xj)i<j<J in C^([0, cx))) with compact support such that 



(p{t, a, x) — ^ ipj {t, x) Xj (a) for (i, a, x) G [0, T) x [0, 00) x . 



(4.12) 



We next use the classical fact that, given G C^([0, T) x [0, 00) x Vt) with compact support, there 
is a sequence of functions {}fn)n which is bounded in Loo{{Q,T) x (0, cxd) x J7) and converges a.e. 
towards </3, each function (y3„ being of the form (|4.12p . The claim (|4.1ip then follows with the help 
of the convergence (|4.6p . 

We next turn to the (strong) compactness of (wq,^ ) ^ and let v denote the solution to 

dtv{t,x) = ig-0{v{t,x)) v{t,x)+Mb^{t,x), (t, x) G (0, T) x 17 , (4.13) 
v{0,x) = v°{x), xen. (4.14) 

Owing to {hi), (/14), and the non-negativity of u and w*^, we have v > and dtv < \\g\\oo v + (3q A. 
Since A belongs to Loo{(0, T)xfl) by (|4.2p and (|4.9p . so does u by the previous differential inequality. 
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It next follows from (hi), (g^l), and that 



1 d 

2 dt 



^2 < 



< 



-a-k llV-WaJIs - Oik 



AVa^ V dx 



2 



M 



M 



bfi 



\{9 - O(Vak) Vak - id - Oiv) V\ \Va^-v\ dx 

2 



-Oik 



Aiiqj. V dx 



M 



M 



+ (1 + ll.9ll=o + C3(T) y - ri|L^(0,C3(T))) \\Va, ^ vf, . (4.15) 
Recalling (|3.27p . we realize that {akAva^)k is weakly relatively compact in L2((0,T) x fl). Since 
{otkVak)k converges to zero in L2((0,T) x fi) by (|4.2p . we thus conclude that {akAva^)k converges 
weakly to zero in L2{{0,T) x fl). Consequently, as v belongs to ioo((0,T) x fi), we have 

lim Uk [ [ Ava^ V dxdt = . (4.16) 

We then infer from (|i?TUl) . (jiTTC)) . and that 

— >v in L2((0,T)xf]) and a.e. in {0,T)xn, (4.17) 

the almost everywhere convergence being obtained after possibly extracting a further subsequence. 

We are now in a position to pass to the limit as a*; ^ in the identity of Lemma 14.21 which 
we first formulate in a different way: ii (p € C^{[0,T) x [0, oo) x Ct) is compactly supported and 
satisfies d,^(p — on [0,T) x [0, oo) x d^, it follows from Lemma l4T2l that 



/ Lp{0,a,x)u°^^{a,x)dadx ^ / / / dtf u^^ dxdadt + 'S2Gn.k{f) 
Jn Jo Jo Jn „^]^ 



with 

Gi,k{^) 

G2MV) 



Jn 



1 



^0 

T 



"fc Jo 

ip{t, a + ak, x) — (p{t, a, x) 
ak 



(4.18) 



93 da ) dxdt , 

- Lp{t, a, x)fiak (a) ) Ua^ {t, a, x) dxdadt , 



I j I Ax^Da^{Ka^)ua^^ dxdadt , 
Jo Jo Jn 

T 



First, the boundedness (|4.2p of {va^)k, the convergence (|4.17p . the continuity of ip, (hi), and the 
Lebesgue dominated convergence theorem allow us to pass to the limit in {Gi^k{^))k and conclude 
that 



lim Gi,k{(p) 

k — >OQ 



Jn 



^{v{t, x)) v{t, x) ip{t,0, x) dxdt . 



(4.19) 



In order to handle {G2,k{^))k, we recall the following consequence of the Dunford-Pettis and Egorov 
theorems, which is implicitly contained in 1 , p. 341] (for a proof see |131 Lem.A.2] for instance). 

Lemma 4.5. Let U he an open hounded subset o/R™, m > 1, and consider two sequences {yk)k G 
Li{U) and {zk)k £ Loc{U) and Junctions y G Li{U) and z G ioo(C^) such that 

Vk ^ y in Li{U) , |zfc(a;)| < K and lim Zk{x) — z{x) a.e. in U 



for some K > Q. Then {ykZk)k converges weakly towards yz in Li{U). 
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We now note that 

/Lp{t,a + ak,x) - ip{t,a,x) \ f \\ a r-^ \ i4- \ ( \ 

hm ip[t, a, xjfia^ (a) = Oaf[t, a, x) - Lp[t, a, x) ^(a) 

fe^oo \ ak / 

for a.e. {t, a, x) G (0, T) x (0, oo) x fl and is bounded in (0, T) x (0, oo) x 17 by (/i4). Due to this 
fact and the weak convergence (j4.6p of {uai^)k in Li we may apply Lemma 14.51 and deduce that 

lim G2 kif) ^ / / / (dafitjajx) — ^p{t,a,x)iJ.{a)) u{t,a,x)dxdadt . (4.20) 
>'^°° ' Jo Jq Jn 



Similarly, {ip Da,^{Aai^))f^ is bounded and converges a.e. in (0, T) x (0, 00) x 17 by virtue of 
and (14. 9|) . so that the same argument applies to establish that 

lim G'3,fc((p) = / / / ^x'pD{A)u dxdadt. (4.21) 
' Jo Jo Jn 

Finally, 3a^ also reads 

= ^(A„j v,.ci(A„j - ^^trr^r^ v.c2(a„j, 

and we infer from (j4.3p and (|4.9p that, after extracting a further subsequence if necessary, we may 
assume that 

V.Ci(A„J-- V,Ci(A) and V,C2(A„ J - V,C2(A) in i2((0, T) x 17) . 
In addition, (^15) and (/ig) imply the local boundedness of D' /([ and -E/C21 so that 

— (A„J^— (A) and ^FTkV ^4((0,r)xl7) 

by (|4.2p . (|4.9p . and (I4.17p . We then conclude from the above two convergence results that 

Jo. -^J:=|^(A)V.Ci(A)-^|^V.C2(A) in L4/3((0, T) x 17) . 

Since ^ 

G4,k{'p) ^ / / Jaj, • A4v,(p,Qfc dxdt 
Jo Jn 

and Mv^-p,ak — * -^v^cp in L4((0,T) x 17) by (|4.1ip . we finally obtain 

lo In {^^^^ ^-^i(^) - ^^r^ V,C2(A)) ■ (^J^ u V.^da) dxdt . (4.22) 



lim G'4,fe(<y3) = 



— >C30 



Due to (|4.19p - (l4.22p . we may pass to the limit as fc ^ 00 in (|4.18p and deduce that u solves (|l.ip 
in the weak sense stated in Definition 12. II □ 



Remark 4.6. It follows from the proof in Section \4.3\ that the strong compactness (j4.17p of {v^) 
and the possibility to pass to the limit in the weak formulation of equation (jl.ip heavily rely on 
the strong compactness (|4.8p of the age averages of{ua)- In turn, this compactness property stems 
from the positivity of D on (0, cxd). Therefore, the interesting case where swarming is only due to 
the drift term diva;(u-E(A, t))Va;A) (corresponding to D = 0) cannot be handled in a straightforward 
way. A similar remark applies to the case where D vanishes on a neighborhood of (e.g. as in 

m)- 



Remark 4.7. It is quite clear from the proof of Theorem \2.2\ that the regularity assumptions 
on the data D, b, X, and the initial data vP , are mainly needed to apply the results in [1] 
for the approximation p.4p - (j3.9p and thus could be weakened. For instance, it would he suf- 
ficient for to be in Loo{^) instead of Wp{Q,). Similarly, one could replace the assumption 
vP e Li((0, 00), VFp(17); A(a)da) by vP G Loo(17; Li((0, 00); A(a)da)). The only modifications to be 
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done in the proof of Theorem \2.2\ would be the construction of suitable approximations {Da), {ha), 
(Aq), {u^a)' '^^^ (^q) ^' ' '^'^'^ ^''j which can be done in a classical way. 
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